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Abstract 

The perturbation series for electroweak vector boson production at small 
transverse momentum is dominated by large double logarithms at each order 
in perturbation theory. An accurate theoretical prediction therefore requires a 
resummation of these logarithms. This can be performed either directly in trans- 
verse momentum space or in impact parameter (Fourier transform) space. While 
both approaches resum the same leading double logarithms, the subleading loga- 
rithms are, in general, treated differently. We comment on two recent approaches 
to this problem, emphasising the particular subleading logarithms resummed in 
each case and the numerical differences in the cross sections which result. 
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A complete theoretical description of W and Z boson production at high-energy 
hadron colliders is necessary for precision Standard Model phenomenology (for exam- 
ple, the measurement of M\y) and for the reliable estimation of backgrounds to new 
physics processes. An important issue is the calculation of the transverse momentum 
{qr) distribution at small qr, which requires the resummation of large 'Sudakov' double 
logarithms. These order dg\r\ 2n ~ 1 {Q 2 / q 2 ^) contributions arise from soft gluon emission 
and lead to a breakdown of fixed-order perturbation theory as qr — > 0. 

Although the resummation of the soft gluon contributions is achieved most natu- 
rally in impact parameter (Fourier transform) space flU , there are certain advantages in 
performing the resummation directly in transverse momentum space 0. For example, 
the matching with fixed-order results at large qr is difficult in the impact parameter 
approach, since this gives oscillatory behaviour for the large qr 'tail' of the transformed 
resummed soft gluon contribution. 

The resummation of the large logarithms directly in qr space has therefore received 
much recent attention. In particular, two complementary approaches have been pro- 
posed in Refs. [fj (FNR) and M (KS). The aim of both these approaches is to improve 
the theoretical description of the small qr cross section by including certain subleading 
(i.e. order a n s \n 2n ~ r {Q 2 / q 2 ^) , r > 2) logarithms in the resummation. However the 
resulting predictions are qualitatively very different at very small qr- The aim of this 
note is to present and discuss the differences between the FNR and KS calculations.^ 
We pinpoint the exact set of subleading logarithms summed in the two cases, and show 
how their inclusion can lead to significant differences in the predictions. 

In order to highlight the differences we will consider both approaches in their most 
simplified version, i.e. we shall restrict ourselves to the parton-level subprocess cross 
section and ignore certain other subleading corrections, as decribed below. We also do 
not consider any non-perturbative {qr smearing) effects. 

The springboard for both approaches is the general expression in impact parameter 
(6) space for the vector boson transverse momentum distribution in the Drell-Yan 
process |] at the quark level: 
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1 We also briefly comment on the related work of g] (EV). 



1 



For the purposes of this simplified analysis we take the coupling as be fixed and retain 
only the leading coefficient i.e. 

A(a s ) = — , B(a g ) = (3) 

7T 

with Cp = 4/3. With these assumptions the b— space expression becomes 

21,2- 
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We recall that the b— space formalism takes fully into account the conservation of 
transverse momentum in multigluon emission, 5(qr + J2i ^n), and is therefore expected 
to provide a better approximation of the q F distribution at small qx than the Double 
Leading Logarithm Approximation (DDLA), in which strong ordering of the gluons' 
kfi is assumed. The DLLA leads to the Sudakov form factor expression || 
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The b— space expression is (mathematically) well defined for all values of q?. In partic- 
ular, in the limit — > it gives a finite positive cross section, in contrast to the DLLA 
result which, because of the vanishing of strong-ordered phase space, yields zero in this 
limit. However the large behaviour of the b— space expression (particularly when 
higher coefficients are taken into account) is not physical - the qx distribution 

oscillates between positive and negative values due to the nature of the Bessel function. 

After performing partial integration the expression (|]) can be written (using the 
KS notation) as 




d6Ji(6)exp --In 2 — j (6) 
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with b = bq T , A = V = $- 

The next step in the KS approach is to expand terms under the integral in (^), 
which gives 

1 da 1 ~ „ (-l/2)«-' ^ (2N-l \ 2N + m 

^ = vL x t^t,£ 2 H » r • (T) 

Here L = ln = ln (jp-^J an d the numbers r m are defined by 

r m = r dyJM^il) (8) 
Jo b 

2 The r m are called b m (oo) in 0. 
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and can be calculated explicitly using the generating function 
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so that e.g. r = 1, t\ = r 2 = 0, r 3 = — |C(3) etc. For large m, the coefficients r m 
behave (— l)' m m!2 m , and the first twenty r m are tabulated in [Q. 

By extracting the Sudakov form factor exp (— f £ 2 ) the expression (0) transforms 



into 
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Naturally, for numerical calculations based on the expression (|ToD it is necessary to in- 
troduce a cut-off value iV max . For example, Fig. [L] shows the contributions to the double 
summation in (10|) which are summed when iV max = 4. Some illustrative numerical 
results based on the KS approach will be presented below. 

Another approach, suggested in (EV), succeeds in developing an analytic ap- 
proximation when a slightly modified set of assumptions is considered. That is to 
say, all the r m coefficients except r are set to be zero and acquires an additional 
contribution — At^{A^) 2 = 2£(3) (A^) 2 . Including the expanded Sudakov factor, see 
Fig. 0, this corresponds to fully summing the first three leading series of logarithms, 
i.e. terms of the form a<V £ 2Ar - 1 -' m ) m = o, 1,2. In the EV approach, the redefinition 
of correctly takes account of the first term of the fourth series, i.e. the 
term. On the other hand, it distorts other terms of this series and terms from more 
subleading series wherever the B^ coefficient appears. 

In the FNR approach 0, one expands the exponent in (|[) 



exp 
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where = In (j^), and retains only the first two terms ('NLL approximation'): 
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Note that keeping only the leading ~ L 2 term in the exponent corresponds to the 
DLLA. A great advantage of the 'NLL approximation' is that the b integral can be 
performed analytically 



1 da 
o"o drj 



[3] 



drj 



exp 




-2AL 



r(i - xl) 
r(i + AL) 



(13) 



is This is 



3 Note that throughout this note the value of the lower limit of integration in 
different from Q where ^ is chosen. Therefore the expression ( |l3| ) differs from the original expression 
in p| by a constant. 
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There is in fact a direct link between the KS and FNR approaches. If instead of 
performing the integration in fll2"D one expands the b— dependent terms in the exponent 
and then performs the integration, the result is 



1 da 
<7 drj 



& rr iL *t^w LN ~ 1[2T « + ^ v (14) 



Clearly this is just the expression fllPP taken at m = 0. Indeed the same result can be 
derived from the resummed expression (|H|) by recalling the definition of the generating 
function and using the relation 

lnr(l + x) = -ln(l + x)+x(l- 7£ ;) + ^(-l) n [C(n)-l]— , |x|<2. (15) 

n=2 71 

The contributions being resummed in both approaches are illustrated schematically in 
Fig. [I] and Fig. 0. The aim of Fig. [I] is to show which terms of the form a^L M from 
the residual sums in ([H]) and (tL4|) are taken into account. Here the FNR approach 
corresponds to having two infinite lines of points (terms) while the KS approach results 
in a finite triangle of terms with size determined by iV max . Terms emerging in the full 
perturbative expansion, i.e. after expanding and multiplying in the Sudakov factor, in 
both approaches are illustrated in Fig. ^. Note that summing over all logarithmic terms 
with a given power of as must result in the perturbative expansion coefficient of the 
same order, up to logarithmic accuracy. Of course a formula with an expanded Sudakov 
factor is valid only when otsL 2 <, 1. The only reason for expanding the Sudakov factor 
here is to determine which terms in the overall perturbation series are actually being 
resummed in (|H]) and ([Hf). It is these latter expressions, which can be regarded as 
the 'master equation' of the two approaches, that we use to obtain numerical results, 
and both approaches remain well-behaved provided a>sL < 1. 

We next present some simple numerical comparisons using both approaches 0. First 
we investigate the dependence of the KS result on the point of truncation iV max . (Note 
that iV max = 1 is just the DLLA approximation.) It is clear from Fig. ^ that for small 
values of r\ (i.e. qr/Q <C 1) the approximating of the b— space result improves with 
increasing iV max . The 'Sudakov dip' at very small 77 is more and more filled in as iV max 
increases, by contributions which are formally subleading in terms of powers of as and 
L. 

On the other hand, the range of applicability of the FNR resummed formula ([13]) is 
seriously restricted. As pointed out in ||, the expression (|13|) suffers from singularities 
at XL = 1, 2, .... (In fact these singularities are poles of order two.) Therefore the first 
pole encountered as rj decreases is at r] cnt = exp ( — x)' It^ = Q ex P (~ 2aJc F )- 
4 The value of the strong coupling constant has been fixed here at as=0.2. 

5 This may appear to be an irrelevantly small value but, as shown in ||, when the running coupling 
constant is used the pole moves significantly towards higher values of qt- 
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Figure 1: Schematic representation of contributions to (0) and to (|10| ). Circles cor- 
respond to the former expression, triangles to the latter one. An empty marker of a 
certain shape means that there exist other contributions in the perturbation series with 
the same power of as and ln(Q 2 /q^) which are not included in an expression coded 
with that shape. The points along the line labelled l M = 2N' represent terms coming 
from the Sudakov factor. 
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Figure 2: Schematic representation of contributions to (|14D with the Sudakov factor 
expanded and (0). Circles correspond to the former expression, triangles to the latter 
one. An empty marker of a certain shape means that there exist other contributions 
in perturbation theory of the same power of as and ln(Q 2 /q^) which are not included 
in an approach coded with that shape. 
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Figure f| shows the resummed FNR result as a function of r]. The pole at ?7 cnt is 
evident (the distribution — > —00 as the singularity is approached from above). The 
resummed result is also compared to the 'truncated' expression (14) for various values 



of the cut-off parameter iV max . This shows the effect of successively adding more and 
more of the subleading 'm = 0' terms along the two infinite lines of Fig. [TJ, starting 
from the Sudakov expression for iV max = 1. Convergence to the (singular) resummed 
FNR result (|13D for large iV max is clearly evident. 

A natural extension of the approach of [[| would be a resummed analytic expression 
also including m = 1 terms in the classification of (|i~4f) . In fact one can systematically 
include the subleading 'NNLL' terms of (|TT|) using the identity 



/A \ 00 1 / A \ J d 2 i 

exp (2LU + £J)] =g - (--) w exp (—XL Li) (16) 

which generates more subleading terms as derivatives of the FNR analytic 'NLL' result. 
In particular, including the m = 0, 1 contributions yields 



1 da 
(T dr] 



m=0A 




2\- 2AL r(i- XL) 
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(17) 



This, however, does not cure the singularity problem but makes it even worse. It turns 
out that if the upper limit of the sum over m increases by 1, the order of the poles 
increases by 2, e.g. the formula flT7| ) has poles of order four at rj cnt . Moreover, as 
this upper limit increases, the region where the approximation of the b— space result 
becomes better, contracts. This is illustrated in Fig. [|, where we show the ratio of the 
numerically calculated expression (|T3|) , (|I~7D and the 'full' b— space result. 

The authors of || argue that the subleading terms in the original expression 
possess a divergent behaviour. These terms have been shown to have factorially grow- 
ing coefficients and originate from the small-6 region of integration. Their presence 
manifests itself in the KS approach which takes some of those subleading terms into 
account. Indeed for large iV max (i.e. more subleading terms), the presence of factorially 
growing subleading coefficients can be seen at large rj, e.g. in Fig. |3|, where the cross 
section can become negative. But this behaviour is not entirely unexpected. Let us 
recall that the b— space formalism was invented to provide a good description in the 
small rj regime where a,sln(l/?7) becomes large. Thus the presence of factorially grow- 
ing terms which manifest themselves for large 7] can be understood as an artefact of 
the b— space method. In fact, in this formalism, the recovery of a credible theoretical 
result in the large rj domain relies on careful matching with fixed-order perturbation 
theory. 

In summary, we have shown that the KS and NFR approaches start from the same 
expression for the cross section in b space, but organise the perturbative expansion in 
qr space in different ways such that different subleading qx logarithms are included. 



7 




10 10 10 10 10 10 10 10 

7] 



Figure 3: The b— space result compared to the expression (|I0|), calculated for various 
values of iV max . Here N max = 1 corresponds to the DLLA approximation. 
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Figure 5: The ratio of the numerically calculated (|13D (m=0 curve) and fllTf) (m= 
curve) to the b— space result. 
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In the KS approach, 'towers' of subleading logarithms fill in the Sudakov (DLLA) dip 
at small q?- In the FNR approach, a particular 'subset' of subleading logarithms is 
resummed to all orders, but the resulting expression has a singularity at q? = qr^i 
below which the cross section is not defined. The FNR result can be obtained in the 
KS approach by including only the appropriate subleading terms. 

The basic question remains as to whether there is a definite value of qx below which 
the perturbative expression for the cross section cannot be calculated. The original 
argument for the 6-space approach was that it allowed a non-zero cross section at qx = 
to be generated by the emission of soft gluons whose transverse momentum vectors 
cancelled, a phase space region clearly outside the strongly-ordered DLLA domain. 
The KS approach is designed to take these (formally subleading) contributions into 
account in a systematic way. The 'price' one pays is a series with factorially growing 
coefficients that drive the behaviour at large q?, but this is in any case outside the 
region of applicability of the whole approach. In contrast, the validity of the FNR 
approach seems to depend on the extent to which the 'LL' and 'NLL' terms as defined 
in Eq. ( |TT| ) do actually give the dominant contribution to the small qx cross section. 
Since we have shown that attempts to systematically include the 'NNLL' contributions 
in this approach lead to even more singular behaviour than the one observed in the 
'NLL' case, this may cast doubt on the validity of the NLL approximation. In any 
case, one incontrovertible conclusion is that this is an interesting and important issue 
that deserves further study. 
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